Measurement-based quantum computation beyond the one-way model 
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We introduce novel schemes for quantum computing based on local measurements on entangled re- 
source states. This work elaborates on the framework established in [Phys. Rev. Lett. 98, 220503 (2007), 
|quant-ph/0609 149") . Our method makes use of tools from many-body physics - matrix product states, finitely 
correlated states or projected entangled pairs states - to show how measurements on entangled states can be 
viewed as processing quantum information. This work hence constitutes an instance where a quantum informa- 
tion problem - how to realize quantum computation ~ was approached using tools from many-body theory and 
not vice versa. We give a more detailed description of the setting, and present a large number of new examples. 
We find novel computational schemes, which differ from the original one-way computer for example in the way 
the randomness of measurement outcomes is handled. Also, schemes are presented where the logical qubits 
are no longer strictly localized on the resource state. Notably, we find a great flexibility in the properties of the 
universal resource states: They may for example exhibit non- vanishing long-range correlation functions or be 
locally arbitrarily close to a pure state. We discuss variants of Kitaev's toric code states as universal resources, 
and contrast this with situations where they can be efficiently classically simulated. This framework opens up a 
way of thinking of tailoring resource states to specific physical systems, such as cold atoms in optical lattices or 
linear optical systems. 



I. INTRODUCTION 

Consider a quantum state of some system consisting of 
many particles. This system could be a collection of cold 
atoms in an optical lattice, or of atoms in cavities, coupled by 
light, or entirely optical systems. Assume that one is capable 
of performing local projective measurements on that system, 
however there is no way to realize a controlled coherent evolu- 
tion. Can one perform universal quantum computing in such a 
setting? Perhaps surprisingly, this is indeed the case: The one- 
way model of Refs. 10, □[ demonstrates that local measure- 
ments on the cluster state - a certain multi-particle entangled 
state on an array of qubits fl] - do possess this computational 
power. The insight gives rise to an appealing view of quantum 
computation: One can in principle abandon the need for any 
unitary control, once the initial state has been prepared. The 
local measurements - a feature that any computing scheme 
would eventually embody - then take the role of preparation 
of the input, the computation proper, and the read-out. This 
is of course a very desirable feature: Quantum computation 
then only amounts to (i) preparing a universal resource state 
and (ii) performing local projective measurements [2-6]. 

But what about other entangled quantum states, different 
from cluster or graph states jsl 0]? Can they form a re- 
source for universal computation? Is it possible to tailor re- 
source states to specific physical systems? For some experi- 
mental implementations - e.g., cold atoms in optical lattices 
lit], atoms in cavities optical systems [11-13], ions in 

traps ifTi]] . or many -body ground states - it may well be that 
preparation of cluster states is unfeasible, costly, or that they 



are particularly fragile to finite temperature or decoherence 
effects. Also, from a fundamental point of view, it is clearly 
interesting to investigate the computational power of many- 
body states - either for the purpose of building measurement- 
based quantum computers or else for deciding which states 
could possibly be classically simulated lfl5i[l6ll . Interestingly, 
very little progress has been made over the last years when 
it comes to going beyond the cluster state as a resource for 
measurement-based quantum computation (MBQC). To our 
knowledge, no single computational model distinct from the 
one-way computer has been developed which would be based 
on local measurements on an algorithm-independent qubit re- 
source state. 

The apparent lack of new schemes for MBQC is all the 
more surprising, given the great advances that have been 
made toward an understanding of the structure of cluster state- 
based computing itself. For example, it has been shown 
that the computational model of the one-way computer and 
teleportation-based approaches to quantum computing ITtIi 
are essentially equivalent ifisi [l9ll . A particularly elegant 
way of realizing this equivalence was discovered in Ref. |20l]: 
They pointed out that the maximally entangled states used 
for the teleportation need not be physical. Instead, the role 
can be taken on by virtual entangled pairs used in a "va- 
lence bond" ll2lll description of the cluster state. This point 
of view is closely related to our approach to be described be- 
low. Further progress includes a clarification of the temporal 
inter-dependence of measurements ||22|] . In Ref. ll23ll a first 
non-cluster (though not universal, but algorithm-dependent) 
resource has been introduced, which includes the natural abil- 
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ity of performing three-qubit gates. Recently, Refs. |]24 
initiated a detailed study of resource states which can be used 
to prepare cluster states (see Section HTAl i. 

In this work, we describe methods for the systematic con- 
struction of new MBQC schemes and resource states. This 
continues a program initiated in Ref. [jTtl in a more detailed 
fashion. We analyze MBQC in terms of "computational ten- 
sor networks", building on a familiar tool from many-body 
physics known by the names of matrix-product states, finitely 
correlated states ll26l l27ll or projected entangled pair states 
MM- 

The problem of finding novel schemes for measurement- 
based computation can be approached from two differ- 
ent points of view. Firstly, one may concentrate on the 
quantum states which provide the computational power of 
measurement-based computing schemes and ask 

1 . What are the properties that render a state a universal 
resource for a measurement-based computing scheme^ ? 

Secondly, putting the emphasize on methods, the central ques- 
tion becomes 

2. How can we systematically construct new schemes for 
measurement-based quantum computation? Is there a 
framework which is flexible enough to allow for the con- 
struction of a variety of different models? 

Both of these intertwined questions will be addressed in this 
work. 



II. MAIN RESULTS 

As our main result, we present a plethora of new universal 
resource states and computational schemes for MBQC. The 
examples have been chosen to demonstrate the flexibility one 
has when constructing models for measurement-based com- 
putation. Indeed, it turns out that many properties one might 
naturally conjecture to be necessary for a state to be a uni- 
versal resource can in fact be relaxed. Needless to say, the 
weaker the requirements are for a many-body state to form a 
resource for quantum computing, the more feasible physical 
implementations of MBQC become. 

Below, we enumerate some specific results concerning the 
properties of resource states. The list pertains to Question 1 
given in the introduction. 

• In the cluster state, every particle is maximally entan- 
gled with the rest of the lattice. Also, the localizable 
entanglement |[29ll is maximal (i.e. one can determin- 
isticaUy prepare an maximally entangled state between 
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' Clearly, the answer to the previous question depends on the definition of 
a universal resource. See Section III Al for a discussion, in particular in 
relation to Ref O. 



FIG. 1: Measurement-based quantum computing as generalization 
of the one-way model as being considered in this work. Initially, 
an entangled resource state is available, different from the cluster 
state, followed by local projective measurements on all individual 
constituents in the regular not necessarily cubic lattice. In all figures, 
dark gray circles denote individual physical systems. 



any two sites, by performing local measurements on the 
remainder). While both properties are essential for the 
original one-way computer, they turn out not to be nec- 
essary for computationally universal resource states. To 
the contrary, we construct universal states which are lo- 
cally arbitrarily pure. 

• For previously known schemes for MBQC, it was es- 
sential that far-apart regions of the state were uncorre- 
cted. This feature allowed one to logically break down 
a measurement-based calculation into small parts cor- 
responding to individual quantum gates. Our frame- 
work does not depend on this restriction and resources 
with non-vanishing correlations between any two sub- 
systems are shown to exist. This property is common 
e.g., in many-body ground-states. 

• Cluster states can be prepared step-wise by means of a 
bi-partite entangling gate (controlled-phase gate). This 
property is important to the original universality proof. 
More generally, one might conjecture that resource 
states must always result from an entangling process 
making use of mutually commuting entangling gates, 
also known as a unitary quantum cellular automaton 
Isoll . Once more, this requirement turns out not to be 
necessary. 

• The cluster states can be used as universal preparators: 
Any quantum state can be distilled out of a sufficiently 
large cluster state by local measurements. Once more, 
this property is essential to the original one-way com- 
puter scheme. However, computationally universal re- 
source states not exhibitingthis properties do exist (the 
reader is referred to Ref. I24{| for an analysis of resource 
states which are required to be preparators; see also the 
discussion in Section HTAl i. More strongly, we construct 
universal resources out of which not even a single two- 
qubit maximally entangled state can be distilled. 

• A genuine qu-trit resource is presented (distinct, of 
course, from a qu-trit version of the cluster state IBlll ). 

We will further see that there is quite some flexibility con- 
cerning the computational model itself (addressing Question 2 
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mentioned in the introduction): 

• The new schemes differ from the one-way model in 
the way the inherent randomness of quantum measure- 
ments is dealt with. 

• We generalize the well-known concept of by-product 
operators to encompass any finite group. E.g. we show 
the existence of computational models, where the by- 
product operators are elements of the entire single-qubit 
Clifford group, or the dihedral group. 

• We explore schemes where each logical qubit is en- 
coded in several neighboring correlation systems (see 
Section|lII]for a definition of the term "correlation sys- 
tem"). 

• One can find ways to construct schemes in which inter- 
actions between logical qubits are controlled by "rout- 
ing" the qubits towards an "interaction zone" or keeping 
them away from it. 

• In many schemes, we adjust the layout of the mea- 
surement pattern dynamically, incorporating informa- 
tion about previous measurement outcomes as we go 
along. In particular, the expected length of a computa- 
tion is random (this constitutes no problem, as the prob- 
ability of exceeding a finite expected length is exponen- 
tially small in the excess). 



A. Universal resource states 

What are the properties from which a universal resource 
state derives its power? After clarifying the terminology, we 
will argue that an answer to this question - desirable as it may 
be - faces formidable obstacles. 

Quantum computation can come in a variety of different 
incarnations, as diverse as e.g., the well-known gate-model 
ll32ll . adiabatic quantum computation Issll or MBQC. All these 
models turn out to be equivalent in that they can simulate each 
other efficiently. 

For measurement-based schemes, the "hardware" con- 
sists of a multi-particle quantum system in an algorithm- 
independent state and a classical computer The input is a 
gate-model description of a quantum computation. In every 
step of the computation, a local measurement is performed on 
the quantum state and the result is fed into the classical com- 
puter Based on the outcomes of previous steps, the computer 
calculates which basis to use for the next measurements and, 
finally, infers the result of the computation from the measure- 
ment outcomes |@|. Having this procedure in mind, we call a 
quantum state a universal resource for MBQC, if a classical 
computer assisted by local measurements on this states can 
efficiently predict the outcome of any quantum computation. 

The reader should be aware that another approach has re- 
cently been described in the literature. The cluster state has 
actually a stronger property than the one just used for the 
definition of universality: it is a universal preparator This 



means that one can prepare any given quantum state on a 
given sub-set of sites of a sufficiently large cluster by means 
of local measurements. Hence, cluster states could in prin- 
ciple be used for information processing tasks which require 
a quantum output. Ref. ll25ll referred to this scenario as CQ- 
universality - i.e. universahty for problems which require a 
classical input but deliver a quantum output. This observation 
is the basis of Ref. ll24ll . where a state is called a universal re- 
source if it possesses the strong property of being a universal 
preparator, or, equivalently, of being CQ-universal. 

Clearly, any efficient universal preparator is also a compu- 
tationally universal resource for MBQC (since one can, in par- 
ticular, prepare the cluster state). But the converse is not true, 
as our results show. Indeed, while it proves possible to come 
up with necessary criteria for a state to be a universal prepara- 
tor ll24ll . we will argue below that the current limited under- 
standing of quantum computers makes it extremely hard to 
specify necessary conditions for computational universality. 

In order to pinpoint the source of the quantum speedup, we 
might try to find schemes where more and more work is done 
by the classical computer, while the employed quantum states 
become "simpler" (e.g., smaller or less entangled). How far 
can we push this program without losing universality? The 
answer is likely to be intractable. Currently, we are not aware 
of a proof that quantum computation is indeed more powerful 
than classical methods. Hence, it can presently not be ex- 
cluded that no assistance from a quantum state is necessary at 
all. 

Observation 1 (Any state may be a universal resource). If 

one is unwilling to assume that there is a separation between 
classical and quantum computation (i.e., BPP ^ BQP), then 
it is impossible to rule out any state as a universal resource. 

It is, however, both common and sensible to assume supe- 
riority of quantum computers and we will from now on do so. 
Observation[T]still serves a purpose: it teaches us that the only 
known way to rule out universality is to invoke this assump- 
tion (this avenue was taken, e.g., in Refs. lfl6[[34ll ). 

Observation 2 (Efficient classical simulation). The only cur- 
rently known method for excluding the possibility that a given 
quantum state forms a universal resource is to show that 
any measurement-based scheme utilizing the state can be effi- 
ciently simulated by a classical computer 

Thus, the situation presents itself as follows: there is a tiny 
set of quantum states for which it is possible to prove that 
any local measurement-based scheme can be efficiently sim- 
ulated. On the other extreme, there is an even tinier set for 
which universality is provable. For the vast majority no as- 
sessment can be made. Furthermore, given the fact that rigor- 
ously establishing the "hardness" of many important problems 
in computer science turned out to be extremely challenging, it 
seems unlikely that this situation will change dramatically in 
the foreseeable future. 

We conclude that a search for necessary conditions for uni- 
versality is likely to remain futile. The converse question, 
however, can be pursued: it is possible to show that many 
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properties that one might naively assume to be present in any 
universal resource are, in fact, unnecessary. 



III. COMPUTATIONAL TENSOR NETWORKS 

The current section is devoted to an in-depth treatment of 
a class of states known respectively as valence-bond states, 
finitely correlated states, matrix product states or projected en- 
tangled pairs states, adapted to our purposes of measurement- 
based quantum computing. This family turns out to be espe- 
cially well-suited for a description of a computing scheme. 

Indeed, any systematic analysis of resources states requires 
a framework for describing quantum states on extended sys- 
tems. We briefly compile a list of desiderata, based on which 
candidate techniques can be assessed. 

• The description should be scalable, so that a class of 
states on systems of arbitrary size can be treated effi- 
ciently. 

• As quantum states which are naturally described in 
terms of one-dimensional topologies have been shown 
to be classically simulable lfl5l[l6ll26ll . the framework 
ought to handle two- or higher dimensional topologies 
naturally. 

• The basic operation in measurement-based computation 
are local measurements. It would be desirable to de- 
scribe the effect of local measurements in a local man- 
ner. Ideally, the class of efficiently describable states 
should be closed under local measurements. 

• The class of describable states should include ele- 
ments which show features that naturally occur in 
ground states of quantum many-body systems, such 
as non-maximal local entropy of entanglement or non- 
vanishing two-point correlations, etc. 

The description of states to be introduced below complies 
with all of these points. 

We will introduce the construction in several steps, starting 
with one-dimensional matrix product states. The new view 
on the processing of information is that the matrices appear- 
ing in the description of resource states are taken literally, as 
operators processing quantum information. 



A. Matrix product states 



• Two D-dimensional vectors \L),\B) representing 

boundary conditions. 

The state vector | ^P) of the matrix product state is then given 
explicitly by ~ 

|*)= (^l^nM...Ai[si]|L) |si,...,s„). (2) 

Si 

From now on we will assume that the matrices are site- 
independent: Ai[j] = A[i], so the MPS is translationally in- 
variant up to the boundary conditions. We take the freedom 
of disregarding normalization whenever this consistently pos- 
sible. 

Let us spend a minute interpreting Eq. (|2|i. Assume we have 
measured the first site in the computational basis and obtained 
the outcome si. One immediately sees that the resulting state 
vector |4''(si)) on the remaining sites is again a MPS, where 
the left-hand side boundary vector now reads 

\L'{si))=A[s,]\L). (3) 

Hence the state of the auxiliary system gets changed according 
to the measurement outcome. So we find that the correlations 
between the state of the first site and the rest of the chain are 
mediated via the auxiliary space, which will thus be referred 
to as correlation space in the sequel. 

In the past, the matrices appearing in the definition of \^) 
have been treated mainly as a collection of variational param- 
eters, used to parametrize ansatz states for ground states of 
spin chains |26]. However - and that is the basic insight un- 
derlying our view on MBQC - Eq. ^ can also be read as an 
operator A[si\ acting on some quantum state \L). We will 
elaborate on this interpretation in Section lTlIBI 

In order to translate Eq. (|2]) to the setting of 2-D lattices, 
we need to cast it into the form of a tensor network. Setting 
Li ~ {i\L) and 

O'l^l*), (4) 

we can write Eq. (|2]) as 

D 

{si,...,Sn\^) = ^ Li^A[si]i„^i^...A[Sn]i„_i,%„Rhr,- 

(5) 

While Eq. (|5]l is awkward enough, the 2-D equivalent is 
completely unintelligible. To cure this problem, we introduce 



A matrix product state (MPS) for a chain of n systems of 
physical dimension d (so d~2 for qubits) is specified by 

• An auxiliary D dimensional vector space (D being 
some parameter, describing the amount of correlation 
between two consecutive blocks of the chain), 

• For each system i a set of _D x D-matrices Ai [j] , j € 
{0...d~l}. 



^ There is a reason why the right-hand-side boundary condition | R) appears 
on the left of Eq. 0. In linear algebra formulas, information usually flows 
from right to left: BA\iIj) means "\ip) is acted on by A, then by B". In 
the graphical notation to be introduce later, it is much more natural to let 
information flow from left to light: 




The order in Eq. ^2) anticipates the graphical notation. 
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a graphical notation-' which enables an intuitive understand- 
ing beyond the 1-D case. In the following, tensors will be 
represented by boxes, indices by edges: 



A[s]i,r 



A[s] 



(6) 
(7) 
(8) 



Needless to say, in the equation above, "Z" is the index leav- 
ing the box on the left-hand-side, "r" the right-hand-side one. 
Connected lines designate contractions of the respective in- 
dices. Eq. (|2]i now reads 



(■si, ■ 



I*) = 



L 








A[Sn] 







A single-index tensor can be interpreted as the expansion co- 
efficients of either a "ket" or a "bra". Sometimes, we will in- 
dicate what interpretation we have in mind by placing arrows 
on the edges: outgoing arrows designating "kets", incoming 
arrows "bras" 



=\l), -J^={r\. 



(9) 



Tensors with two indices Ai^r can naturally be interpreted as 
operators. In the graphical notation we often want to think of 
information flowing from the left to the right, in which case 
^ ^ r Ai,rV)j.{l\i would be denoted as 



A 



= A, 



(10) 



i.e. with the l.h.s. index being associated with a "bra" and the 
rh.s one with a "ket". The following relations exemplify the 
definition: 



{R\L) = \l\-{r\, 
A\L) = 



AB 



tr(AB) = 



B 



A 



B 



A 



(11) 



(12) 



(13) 



(14) 



The formula for the expansion coefficients of a matrix product 
state finally becomes 



L 




A[s,] 




A[Sn] 













This formula suggest a more "dynamic" interpretation of 
MPS: the l.h.s. boundary conditions \L) specify an initial state 
of the correlation system, which is acted on by the matrices of 
the MPS representation. The next paragraph is going to elab- 
orate on this point. 



B. Quantum computing in correlation systems 

We return to the discussion of the properties of matrix prod- 
uct states. Above, it has been shown how to compute the over- 
lap of I'l') with an element of the computational basis (c.f. Eq. 
(|5]l). The next step is to generalize this to any local projec- 
tion operator. Indeed, if | 
abbreviate 



is a general state vector in C-, we 



One then easily derives the following, central formula 



(15) 



(0.1)1^) 



L 












R 



(16) 

Now suppose we measure local observables on \^) and ob- 
tain results corresponding to the eigenvector at the i-th 
site. Eq. ( fT6] l allows us to re-interpret this process as follows. 
Initially, the /^-dimensional correlation system is prepared in 
the state \L). The result \(f)i) at the first site induces the evo- 
lution 



\L) 



A[0i]|L). 



(17) 



From this point of view, a sequence of measurements on \^) 
is tantamount to a processing of the correlation system's state 
by the operations An appeahng perspective on MBC 

suggests itself: 

Observation 3 (Role of correlation space). Measurement- 
based computing takes place in correlation space. The gates 
acting on the correlation systems are determined by local 
measurements. Intuitively, "quantum correlations" are the 
source of a resource's computational potency. The strength 
of this framework lies in the fact that it assigns a concrete 
mathematical object to these correlations. 

Indeed, it will turn out that MBQC can be understood com- 
pletely using this interpretation. 



C. Example: The 1-D cluster state 

To illustrate the abstract definitions made above, we will 
discuss the linear cluster state vector \Cln) in this section. It 
is both one of the simplest and certainly the most important 
MPS in the context of MBQC. 

What is the tensor network representation of |C/„)? Recall 
that the cluster state can be generated by preparing n sites in 
the state vector |+) := |0) + |1) and subsequently applying 
the controlled- Z operation 



CZ= |0,0)(0,0| + |0,1)(0,1| + |1,0)(1,0| 



ll,l)(l,ll 

(18) 



^ These graphical formulae are compatible with various similar systems in- 
troduced before (36ll . 



^ Of course, for general measurement bases, A[(^i] is not going to be unitary. 
Choosing the bases in such a way as to ensure unitarity is an essential part 
of the design of a computational scheme for a given resource. 
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between any two nearest neighbors. Effectively, CZ intro- 
duces a TT -phase whenever two consecutive systems are in the 
|l)-state. Hence its expansion coefficients in the computa- 
tional basis are given by 



■ , Sn\Clr, 



1/2, 



(19) 



where p denotes the number of sites i such that Si = Si+i = 1. 

This observation makes it simple to derive the tensors of 
the MPS representation. We need a D ~ 2-dimensional cor- 
relation system, which - loosely speaking - will convey the 
information about the state Si of the i-th site to site i + 1. 
Define the matrices A [0/1] by 



A[0] 



(20) 
(21) 



The intuition behind this choice is as follows. By the elemen- 
tary relations 

(+io) = (+ii) = (-10) = 2-l/^ (-11) = -2-l/^ 

(22) 

the contraction in the middle of 



A[si] — A[s2] 



(23) 



will yield a sign of "— i" exactly if si = S2 = 1. Indeed, 
setting the boundary vectors to \L) = |0), \R) = |+) one 
checks easily that 



{R\A[s„] . . . A[si]\L) 



-n/2 



(-If, 



(24) 



which is exactly the value required by Eq. iT% . 

Below, we will interpret the correlation system of a 1-D 
chain as a single logical quantum system. For this interpre- 
tation to be viable, we must check that the following basic 
operations can be performed deterministically by local mea- 
surements: i) prepare the correlation system in a known initial 
state, ii) transport that state along the chain (possibly subject 
to known unitary transformations) and iii) read out the final 
state. 

To set the state of the correlation system to a definitive 
value, we measure some site - say the i-th - in the Z- 
eigenbasis. Throughout this work, we will choose the notation 
X, Y, and Z for the Pauli operators. Denote the measure- 
ment outcome by z e {0, 1}. In case of z = 0, Eq. ( |20] | tells 
us that the state of the correlation system to the right of the 
i-th site will be |+) (up to an unimportant phase). Likewise, a 
z = 1 outcome prepares the correlation system in | — ) , accord- 
ing to Eq. ( I2TI 1. It follows that we can use Z-measurements 
for preparation. How to cope with the intrinsic randomness of 
quantum measurements will concern us later. 

Secondly, consider the operators 



A[+] ^ = 2-i/2( ^ A[Q] 



A[l] 



« \+){0\ + \-){l\ = H, 



cx HZ, 



' ) 
(25) 

(26) 



where H is the Hadamard-gate. We see immediately that mea- 
surements in the X-eigenbasis give rise to a unitary evolution 
on the correlation space. Similarly, one can show that one 
can generate arbitrary local unitaries by appropriate measure- 
ments in the Y-Z plane. 

Below, we will frequently be confronted with a situation 
like the one presented in Eqs. (|25|26t , where the correlation 
system evolves in one of two possibilities, dependent on the 
outcome of a measurement. It will be convenient to introduce 
a compact notation that encompasses both cases in a single 
equation. So Eqs. ( 125126b will be represented as 



A[X] 



HZ"". 



(27) 



Here x = corresponds to the outcome |+) in an X- 
measurement, whereas x = \ corresponds to the outcome 
|— ). In general, a physical observable given as an argument 
to a tensor corresponds to a measurement in the observable's 
eigenbasis. The measurement outcome is assigned to a suit- 
able variable as in the above example. 

Lastly, we must show how to physically read out the state 
of the purely logical correlation system. It turns out that mea- 
suring the i + 1-th physical system in the Z-eigenbasis cor- 
responds to a Z-measurement of the state of the correlation 
system just after site i. Indeed, suppose we have measured the 
first i systems and obtained results corresponding to the local 
projection operator \4>i) ® ■ ■ ■ ® \(t)i) ■ Further assume that as 
a result of these measurements the correlation system is in the 
state |0): 



L 











|0). (28) 



Using Eq. ( I2TI 1 we have that 



L — 










A\l] 



(29) 



« |+)(1|0)=0. 



But then it follows from Eq. ( fT6b that the probability of ob- 
taining the result 1 for a Z-measurement on site i + 1 is equal 
to zero. In other words: if the correlation system is in the state 
|0) after the i-th site, then the i + 1-th physical site must also 
be in the state |0). An analogous argument for the |l)-case 
completes the description of the read-out scheme. 



D. 2-D lattices 

The graphical notation greatly facilitates the passage to 2- 
D lattices. Here, the tensors A[s\ have four indices A[s\i^r,u,d, 
which will be contracted with the indices of the left, right, up- 
per and lower neighboring tensors respectively. After choos- 
ing a set of boundary conditions |_L), |i?), \U), \D) £ C^, the 
expansion coefficients of the state vector j^*) are computed as 
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illustrated in the following example on a 2 x 2-lattice: 



in Section lTlICI The resulting tensors are: 



U 



U 



(Sl,l,-..,S2,2|«') 



R 



A[s 



1,2 



A[s2,: 



(30) 



R 



D 



D 



In the 1-D case, we thought of the quantum information 
as moving along a single correlation system from the left to 
the right. For higher-dimensional lattices, a greater deal of 
flexibility proves to be expedient. For example, sometimes 
it will be natural to interpret the tensor Ai,r,u.d as specifying 
the matrix elements of an operator A mapping the left and the 
lower correlation systems to the right and the upper ones: 



{r\(g,{u\A\l) (g,\d), A 



Often, on the other hand, the interpretation 



^,d = {r\A\l) (g)\u) (g)\d), A 



(31) 



A 



(32) 



or yet another one is to be preferred. 

We have seen in Section lTlIBI that the correlation system of 
a one-dimensional matrix product state can naturally be inter- 
preted as a single quantum system subject to a time evolution 
induced by local measurements. It would be desirable to carry 
this intuition over to the 2-D case. Indeed, most of the exam- 
ples to be discussed below are all similar in relying on the 
same basic scenario: some horizontal lines in the lattice are 
interpreted as effectively one-dimensional systems, in which 
the logical qubits travel from the left to the right. The vertical 
dimension is used to either couple the logical systems or iso- 
late them from each other (see Fig.[T]i- The reader should re- 
call that this setting is very similar to the original cluster state 
based-techniques. Clearly, it would be interesting to devise 
schemes not working in this way and the example presented 
in Section lTV B 2l takes a first step in this direction. 



E. Example: the 2-D cluster state 

Once again the cluster state serves as an example. One can 
work out the tensor network representation of the 2-D cluster 
state vector \ Clnxn) in the same way utilized for the 1-D case 



Jl 



A[0] 

"T" 
A. 



A[l] 
"T" 



\L) = \D) = \+), \R)^\U) = \1). 



(33) 



(34) 



(35) 



An important property of Eqs. ( l33l [34b is that the tensors 
A [0/1] factor One could graphically represent this fact by 
writing 



where 



-A[Q]- = ^ {+}- 
^ =|0>,^- =|+). 



(36) 



(37) 



In other words: the tensors ^[0/1] effectively de-couple their 
respective indices. Based on this fact, we will see momen- 
tarily how Z-measurements can be used to stop information 
from flowing through the lattice. 

Indeed, suppose three vertically adjacent sites are mea- 
sured, from top to bottom, respectively in the Z, X and Z- 
eigenbasis: 



Jl 



A[Zu 



A[X] 



(38) 



A[Za 



Denote the measurement results by z„,a;, G {0,1}. As 
before, these numbers correspond to 0„ = for |0) and z„ = 
1 for |1), as well as a; = for |+) and x = \ for |— ). In fact, 
we are mainly interested in the indices of the middle tensor, as 
they will be the ones which carry the logical information. To 
this end Eq. ( [36] l is of use, as it says that the upper and lower 
tensors factor and hence it makes sense to dis-regard all of 
their indices which do not influence the middle part. It hence 
suffices to consider 



A[Zu] 



A[X] 



(39) 



A[Zd] 
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As a first step, we calculate 



1 



|0J 

r+1 



A[0] 



^ ^ =2-i|+)(0|, 



having used Eq. ( 1361 ) and the basic fact 

(0|+) =2-1/2. (4Q^ 

A similar calculation where A[0] is substituted by A[l] yields 
2-i|-)(l|. Hence, for A [+] cx A[0] + we have 



A[+] ^ (X |+)(0| + |-)(1| =i7. (41) 



Similarly, 



cx iJZ. 



(42) 



+ 



After these preparations it is simple to conclude that 



A\Z„, 



A[X] 



(43) 



A[Zd 



This finding tells us how to transport quantum information 
along horizontal lines through the lattice. Namely by measur- 
ing the line in the X-eigenbasis to cause the information to 
flow from the left to the right and measuring vertically adja- 
cent sites in the Z-eigenbasis to shield the information from 
the rest of the lattice. 

Eq. (03]) should be compared with Eqs. ( 125126b . So up to 
possible corrections of the form the procedure out- 

lined above enables us to effectively prepare a 1-D cluster 
state within the 2-D lattice. 



It is the purpose of this section to discuss a number of exam- 
ples of new resources. As before, important features will be 
highlighted as "observations". 



A. AKLT-type states 

1. 1-D structures 

Our first example is inspired by the AKLT state ll2lll . which 
is well-known in the context of condensed matter physics. The 
AKLT model is a 1-D, spin-1, nearest neighbor, frustration 
free, gapped Hamiltonian. Its unique ground state is a ma- 
trix product state with D ~ 2 and indeed, the AKLT model 
motivated the first studies of such states l.21,.,26il . The defining 
matrices of the MPS description are: 



^[0] 



^[2] 



= Z, (44) 
= 2-i/2|0),(l|i, (45) 
= 2-i/2|l),(0|i (46) 



We will choose the boundary conditions to be |L) = \R) = 
|0). As a matter of fact, we will not work directly with the 
AKLT state, but with a small variation, for which it turns out 
to be more straight-forward to construct a scheme for MBQC. 
In this modification, the matrix ^[0] is given by the Hadamard 
gate, instead of the Pauli Z operator; 



A[0] 



= H. 



(47) 



This state shares all the defining properties of the original: it 
is the unique ground-state of a spin-1 nearest neighbor frustra- 
tion free gapped Hamiltonian (see Appendix l VIII Bl ). Against 
the background of our program, the obvious question to ask is 
whether these matrices can be used to implement any evolu- 
tion on the correlation space. 

To show that this is indeed the case, let us first analyze 
a measurement in the {|0), |+), |— )}-basis, where |±) := 
2-i/2(|l^ ± |2)). In a mild abuse of notation, we will 
hence write |±) for state vectors in the subspace spanned 
by {|1), |2)} instead of {|0), |1)}. From Eqs. ( I44B7] ) one 
finds that depending on the measurement outcome, the op- 
eration realized on the correlation space will be one of iJ, X 
or ZX ^ lY . At this point, we have to turn to an impor- 
tant issue: how to compensate for the randomness of quantum 
measurement outcomes. 



rV. NOVEL RESOURCE STATES 

Up to this point, we have reformulated the computational 
model of the one-way computer in the language of computa- 
tional tensor networks. This picture of one-way computation 
is educational in its own right. However, to convincingly ar- 
gue that the framework is rich enough to allow for quite dif- 
ferent models, we have to explicitly construct novel schemes. 



2. Compensating the randomness 

Assume for now that we intended to just transport the infor- 
mation faithfully from left to right. In this case, we consider 
the operator 

Bl -.^ H,X, 01 ZX (48) 
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as an unwanted by-product of the scheme. The one-way 
computer based on cluster states has the remarkable prop- 
erty that the by-products can be dealt with by adjusting 
the measurement-bases depending on the previous outcomes, 
without changing the general "layout" (in the sense of Fig. 
[T]) of the computation yD. For more general models, as the 
ones considered in this work, such a simple solution seems not 
available. Fortunately, we can employ a "trial-until-success" 
strategy, which proves remarkably general. 

The key points to notice are that i) the three possible out- 
comes H, X and Z generate a finite group B and ii) the prob- 
ability for each outcome is equal to 1 /3, independent of the 
state of the correlation system. We will refer to B as the 
model's by-product group. Now suppose we measure m ad- 
jacent sites in the {|0), |+), |— )}-basis. The resulting overall 
by-product operator B = B„iB,n_i . . . Bi will be a product 
of m generators H, X, ZX. So by repeatedly transporting the 
state of the correlation system to the right, the by-products are 
subject to a random walk on B. Because B is finite, every ele- 
ment will occur after a finite expected number of steps (as one 
can easily prove). 

The group structure opens up a way of dealing with the 
randomness. Indeed, assume that initially the state vector of 
the correlation system is given by B\ip), for some unwanted 
B G B. Transferring the state along the chain will introduce 
the additional by-product operator B~^ after some finite ex- 
pected number of steps, leaving us with 



B~'B\^j) = IV), 



(49) 



as desired. The technique outlined here proves to be extremely 
general and we will encounter it in further examples presented 
below. 

Observation 4 (Compensating randomness). Possible sets of 
by-product operators are not limited to the Pauli group. A 
way of compensating randomness for other finite by-product 
operator groups is to adopt a "trial-until-success strategy", 
which gives rise to a random length of the computation. This 
length is in each case shown to be bounded on average by a 
constant in the system size. 



3. All single-qubit gates 

By the preceding paragraphs, we can implement any el- 
ement of B on the correlation space. We next address the 
problem of realizing a phase gate S'(0) := diag(l, e**^) for 
some (j> £ M.. To this end, consider a measurement on the 
{|0),2-i/2(|x^ ± e''^|2))}-basis. There are three cases 



• The outcome corresponds to |1) 



1 2). In this case. 



we get 5(0) on the correlation space and are hence 
done. 

The outcome corresponds to |1) — e*'^|2). We get 
ZS{4>), which is the desired operation, up to an element 
of the by-product group, which we can rid ourselves of 
as described above. 



i: 



(OiOi+r) - {Oi) {0,+r) OC e 



• •- • • • 

I I I I I 



FIG. 2: A universal resource deriving from the AKLT-model. 

• Lastly, in case of |0), we implement H on the correla- 
tion space. As H €i B, we can "undo" it and then re-try 
to implement the phase gate. 

Hence, we can implement any element of B as well as S{(f>) 
on the correlation space. This implies that HS{4))H is also 
realizable and therefore any single-qubit unitary, as SU (2) is 
generated by operations of the form S'((/)) and HS{4))H. 

The state of the correlation system can be prepared by mea- 
suring in the computational basis. In case one obtains a re- 
sult of "1" or "2", the state of the correlation system will be 
|0) or |1) respectively, irrespective of its previous state. A 
"0"-outcome will not leave the correlation system in a defi- 
nite state. However, after a finite expected number of steps, 
a measurement will give a non-"0"-result. Lastly, a read-out 
scheme can be realized similarly (c.f. Section lTlI Cl l. 

Observation 5 (Ground states). Ground states of one- 
dimensional gapped nearest-neighbor Hamiltonians may 
serve as resources for transport and arbitrary rotations. 



4. 2-D structures 

Several horizontal 1-D AKLT-type states can be coupled to 
become a universal 2-D resource. The coupling can be facili- 
tated by performing a controlled-Z operation, embedded into 
the three-dimensional spin-1 space, between vertically adja- 
cent nearest neighbors. More specifically, we will use the op- 
eration exp{i7r|2) (2| (g) |2)(2|}, which introduces a 7r-phase 
between two systems exactly if both are in the state |2). The 
tensor network representation of this resource is given by 



^[0] 



m 



2-'/'|0).(l| 



(50) 



Id, (51) 



A[2] 



2-'/'|l).(0|,®h)Jl|„ (52) 
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as one can check in analogy to Sec. lIIIE] Here, 

Hi^r :=|+).(0|; + |-),,(1|;. 



(53) 



To verify that the resulting 2-D state constitutes a universal 
resource, we need to check that a) one can isolate the corre- 
lation system of a horizontal line from the rest of the lattice, 
so that it may be interpreted as a logical qubit and b) one can 
couple these logical qubits to perform an entangling gate. 

The first step works in complete analogy to Section Ull El 
see Fig.|2l Indeed, one simply confirms that 



A[Zu 



A[s] 



± 



(54) 



where s S {0,1,2} and Z^/i denotes a measurement in 
the {|0), |2)}-basis. So measuring the vertically adjacent 
nodes in the computational basis gives us back the 1-D state, 
up to a possible sign. 

A controlled-Z gate can be realized in five steps: 







A[X] — A[X] — A[X] — A[X] — A[X 



A[Z] — A[Z] 



A\Y] 



A[Z] — A[Z] 



A[X] — A[X] — A[X] — A[X] — A[X 



(55) 



The Pauli matrices X, Y, Z are understood as being embedded 
into the {|1), |2)}-subspace. So, e.g., X denotes a measure- 
ment in the {\0) ,2^^/\\l) ± |2))}-basis. When operating the 
gate, we first measure all sites of the upper and lower lines in 
the X-eigenbasis. In case the result for the sites at position 
"0" (refer to labeling above) is different from |+), the gate 
failed. In that case all sites on the middle line are measured 
in the computational basis and we restart the procedure five 
steps to the right^. Otherwise, the systems labeled by a Z are 
measured. We accept the outcome only if we obtained |1) on 
sites ±2 and |0) on sites ±1 - should a different result occur, 
the gate is once again considered a failure and we proceed as 
above. Lastly, the Y measurement on the central site is per- 
formed. In case of a result corresponding to |0), it is easy to 
see that no interaction between the upper and the lower part 
takes place, so this is the last possibility for the gate to fail. 



Let us assume now that the desired measurement outcomes 
were realized. At site —2 on the middle line, we obtained 



A[l] 



(56) 



which prepares the correlation system of the middle line in 
|0). At site —1, in turn, a Hadamard gate has been realized, 
which causes the output of site —1 to be H\0) = |+). The 
situation is similar on the rh.s., so that the above network at 
site can be re-written as 



A[- 



+ 




A\Y] 




+ 



(57) 



A[- 



We will now analyze the tensor network in Eq. (|57] | step by 
step. For proving its functionality, there is no loss of general- 
ity in restricting attention to the situation where the correlation 
system of the lower line is initially in state |c), for c G {0, 1}. 
We compute for the lower part of the tensor network 



JL 



A[+] 



X\c)Z'' 



(58) 



Further, plugging the output Z'^'\+) of the lower stage into the 
middle part, we find 



JL 



A[Y] 



oc Z^+^'dO) (59) 



Z^\ 



where y G 0, 1 reflects the outcome of the F-measurement on 
the central site: ?/ = in case of |1) + i|2) and ?/ = 1 for 
|1) -i|2). Lastly, 



sz^+yx. 



z^+n|o) + z|i)) 



(60) 



In summary, the evolution afforded on the upper line is 
HSZy^"^, equivalent to Z"^ up to by-products. This completes 
the proof of universality. 

For completeness, note that we never need the by-products 
to vanish for all logical qubits of the full computation simul- 
taneously. Hence the expected number of steps for the real- 
ization of one- or two-qubit gates is a constant in the number 
of total logical qubits. 



^ We have chosen this approach in order to avoid an awkward discussion 
of how to handle phases introduced by "wrong" measurement outcomes. 
We are providing proofs of principle for universality here and will accept a 
(possibly daunting) linear overhead in the expected number of steps, if this 
simplifies the discussion. Substantial improvements to these schemes are, 
of course, possible. 



B. Toric code states 

In the following, we present two MBQC resource states 
which are motivated by Kitaev's toric code states IJs!]. This 
contrasts with a result in Ref. Isill that MBQC on the planar 
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toric code state itself can be simulated efficiently classically. 
Different from the other schemes presented, the natural gate 
in these schemes is a two-qubit interaction, whereas local op- 
erations have to be implemented indirectly. Also, individual 
qubits are decoupled not by erasing sites but by switching off 
the coupling between them. 

Toric code states are states with non-trivial topological 
properties and have been introduced in the context of quan- 
tum error correction. They have a particularly simple repre- 
sentation in terms of PEPS ll39ll or CTNs iQl] on two centered 
square lattices. 



Kh 




Kh 




Kh 



Kv 




Kv 




Kv 



Kh 




Kh 




Kh 



Kv 



Kv Kv 



(61) 



Thus, the tensors in Kitaev's toric code state have a two- 
qubit operation as their natural gate in correlation space, rather 
than a single-qubii gate. In MBQC schemes which base on 
these projectors, two-qubit gates are easy to realize, whereas 
in order to get one-qubit gates, tricks have to be used. In the 
first example, we obtain single-qubit operations by introduc- 
ing ancillae: a ZZ controlled phase between a logical qubit 
and an ancilla in a computational basis state yields a local 
Z rotation on the logical qubit. In the second example, we 
use a different approach: we encode each logical qubit in two 
qubits in correlation space. Using this nonlocal encoding, we 
obtain an easy implementation of both one- and two-qubit op- 
erations; furthermore, the scheme allows for an arbitrary par- 
allelization of the two-qubit interactions. 

Observation 6 (Logical qubits in several correlation systems). 
There is no need to have a one-one correspondance between 
logical qubits and a single correlation system. 



1. Toric codes: first scheme 



where 



Our first scheme consists of the modified tensor 



Kh[s] 



(62) 



Kh[s] 



Kh[s] 



ZH 



(66) 



and 



Kv\s\ 



Z" 



(63) 



ZHZ'' 



i.e., Kh and Kv are identical up to a rotation by 90 degrees. 

Let us first see how Kh acts on two qubits in correla- 
tion space coming from the left. The most basic operation 
is a measurement in the computational basis, which simply 
transports both qubits to the right (up to a correlated Z by- 
product operator). Generalizing this to measurements in the 
Y-Z plane, we find that 



KhW\ 



ZZ{<\>) 



(64) 



where is the angle with the Z axis, and 



/I 



V 



oi0 



(65) 



(Note that this gate is locally equivalent to the CNOT gate for 
(/) = ±7r/2.) 



[with \fZ — diag(l, i)], arranged as in (jSTJ where both Kh 
and Kv are replaced by Kh- The extra H serves the same 
purpose as in other schemes: it allows to leave the subspace of 
diagonal operations and thus to implement X rotations. The 
need for the \fZ will become clear later; it is connected to the 
fact that 



CNOT = {\®H) {VZ ® VZ) ZZ{-tt/2) {1(g) H) . (67) 

In the following, we show how this state can be used for 
MBQC. The qubits run from left to right in correlation space 
in zig-zag lines in Eq. (1611 1; for the illustration in Fig. [3] we 
have straightened these lines, and marked the measurement- 
induced ZZ interactions coming from the Kh [s] in ( |66] l by 
ellipses. (The difference between filled and non-filled ellipses 
will be explained later) The \/ZH operations of ( |66] | do not 
depend on the measurement and are thus hard-wired; note that 
the order is reversed as we are considering H and \/Z as two 
independent operations in the circuit. 

Let us first impose that all qubits are initialized to |0); this 
corresponds to a left boundary condition |0) in correlation 
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a) 



(01 



- 1 h4z I 



(0| V\ // yn 47. 



H^Z 



H^Z 



h4z 



c) 



(01 



\H\\^Z[ 



H^ZI 



h4z[ 



•H^JZ 



" "^"^ancillae 
'^/- log. qubil 



(OH 



ZZ(0) 



log. qubit ^ H ^ — — H — 

zz{4>) 

(o| ^ \ 



]H\^Z 




logical interaction 



iVZ 



iifVZ 



-ffiVZ s 



!//VZ 



■ ancilla 

■ log. qubit 



d) 



PC 



(01- 



Vz 



ZZ(0) 



e) 



— 



ZZ(0) 



H 



H — Vz — 



Z — 

Z — 



FIG. 3: Implementation of single-qubit and two-qubit operations in the first toric code model, a) The measurement pattern for single-qubit 
operations and b) the corresponding circuit, c) Pattern for a two-qubit gate between logical qubits, d) the corresponding circuit and e) the 
circuit after some simplifications. 



space. We will discuss later how to initialize the scheme. Ev- 
ery second qubit is an ancilla which will be used to implement 
one-qubit operations. We first discuss the case of no Pauli er- 
rors, and show later how those can be dealt with. 

The implementation of single-qubit operations is illustrated 
in Fig. |3^. There, each ellipse denotes a possible ZZ inter- 
action. In particular, empty ellipses denote interactions which 
are switched off (i.e. measured in the Z basis), while filled el- 
lipses denote sites where one can measure in the Y-Z plane to 
implement a ZZ gate. If all interactions are switched off, all 
qubits are transported to the right, subject to the transforma- 
tion \^H. As {sfZH)^ = 1, the ancillae are in the compu- 
tational basis in every third step: These regions are hashed in 
Fig. [3^. In these regions, a ZZ{(t>) between ancilla and logical 
qubit (corresponding to the filled ellipses in the figure) results 
in a single-qubit Z rotation on the latter Thus, in each block 
of length three as the one shown in Fig. [3^, the transformation 

y/ZHVZHS{tp)\^HS{(j)) = HS{ip)HS{^) (68) 

is implemented [where S{4i) = diag(l, e**^)], which allows 
for arbitrary one-qubit operations. In Fig. [3j', the corre- 
sponding circuit is shown, which has been simplified using 
H^/ZHVZ = sfX^fZ = {^fzy^H, and that diagonal ma- 
trices commute. 

Although the scheme has a natural two-qubit interaction, 
implementing an interaction between two adjacent logical 
qubits is complicated by the ancilla which is located inbe- 
tween. In order to obtain a coupling, we first swap the logical 
qubit with the ancilla, then couple it to the now adjacent logi- 
cal neighbor, and finally swap it back. This is implemented by 
the measurement pattern shown in Fig. Again, empty el- 
lipses correspond to switched off interactions, while the filled 



ellipses all implement ZZ{—n/2) gates, each of which to- 
gether with two \/Z and two Hadamards as grouped in the 
figure gives a CNOT gate, cf. Eq. ( l67b . This measurement pat- 
tern corresponds to the circuit shown in Fig. [3]l, where we 
have replaced each pair of CNOTs by a CNOT and a SWAP. 
By merging each CNOT with the two adjacent Hadamards, we 
effectively obtain 

CZ = |0,0)(0,0| + |0, 1)(0,1| + |1,0)(1,0| - |1, 1)(1,1| 

(69) 

gates. We thus remain with only diagonal gates on the two 
lower qubits (except for the SWAP), i.e. the gates all commute 
and the circuit can thus be simplified to the one shown on 
in Fig. [3^, proving that the sequence effectively implements 
a two-qubit interaction between the logical qubits. Note that 
the length of the complete sequence is compatible with the 
three-periodicity of the basis of the ancillae. 

Pauli errors in this scheme can be dealt with as usual; H 
and Vz are both in the Clifford group, i.e., Paulis can be 
commuted through, and ZZ commutes with Z errors, while 
(1 ® X)ZZ{(j)) = ZZ{-(t)){l ® X). 

Finally, we show how to read out the logical qubits. It holds 
that 
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i.e., a measurement in the X basis returns the parity of the an- 
cilla and the logical qubit. If this is done when the ancilla is in 
a computational basis state, one effectively measures the logi- 
cal qubit in the computational basis. Note that both the ancilla 
and the logical qubit are in a well-defined state afterwards and 
can thus be reused. 

Let us now turn towards the initialization procedure. In 
contrast to the previous MBQC schemes, the read-out cannot 
be used for initialization. The reason is that the read-out only 
works if the ancilla qubit is initially in a computational basis 
state; otherwise, it just projects onto the subspace spanned by 
{|0,0),|l,l)}orby{|0,l),|l,0)}. 

In the following, we demonstrate that it is still possible 
to initialize this scheme by taking a different perspective on 
how it encodes logical qubits. Therefore, we group each log- 
ical qubit with the ancilla above (e.g., the first two qubits in 
Fig. [3^), and encode the new logical qubit in their parity - 
note that this is what is really measured in the read-out. The 
following calculations are most conveniently carried out in a 
Bell basis where each state is described as |s)^ 1^)^, where the s 
qubit stores the sign of the Bell state and the I qubit the parity 
and thus encodes our logical qubit, i.e. 



k>jO)z - |0,0) + (-1)^11,1) 
|s)Jl), - |0,1) + (-1)11,0) 



(72) 
(73) 



The circuit transforming between the above encoding and the 
qubits in correlation space is 



\s)— H 



ancilla 



^ — logical qubit 



(74) 



commute, and {CZY = 1- As we see from the simplified 
circuit, we obtain an X rotation on the upper logical qubit, 
but with the rotation direction determined by the state of the 
|s)^ qubit below: While |+)^ results in a rotation Rx{4>), the 
state |— )^ gives 

ZR,{(j})Z ^R,{-4,) . 

Similarly, the circuit for the coupling of two logical qubits can 
be simplified as in Fig. 2};: again, the coupling on the logical 
qubits is Cpl((/)) := [H ® Z)ZZ{(f)){H (g) 1) or 

{H (g) ZX)ZZ{(t>){H (g)X) = {Z(g) l)Cpl(0)(Z (g) 1) , 

depending on whether the second s qubit is |+)^, or |— )^. 

Therefore, the error introduced by the unknown state of 
each s qubit results in a Z correction around each operation on 
the logical qubit above (note that we can assume this also for 
Z rotations as they commute with the Z correction). Although 
the error itself is unknown and different for each logical qubit, 
it is consistent within each qubit, as it is always determined by 
the same ancilla. Thus, two subsequent Z errors cancel out, 
and one remains only with one Z correction on the logical 
qubit at the beginning and one at the end of the sequence. The 
former has no effect since the initial state is |0)j, while the 
latter has no effect either since the encoded logical qubit is 
finally measured in the computational basis. Thus, the output 
statistics for the circuit is independent of the initial state |±)g 
of the phase qubits, and one can equally well start from their 
mixture 1^ which completes the argument. 



Using this decoding, it is straightforward to investigate 
what happens in the various steps of the MBQC scheme. 
Firstly, one can easily check that by measuring two consec- 
utive couplings of the qubit pair in the X basis, one prepares 
them in a maximally entangled state |0, 0) + 1 1, 1) up to Pauli 
errors, corresponding to |0)j,|0); in the encoded system. By 
pretending a Pauli Z error on one of the qubits with p = 1/2, 
we effectively face the mixture |0,0)(0,0| + |1, 1)(1, 1|, cor- 
responding to Is (Xi |0)(0|;. 

Since the transformation ( 1741 1 is in the Clifford group, Pauli 
errors remain Pauli errors in the encoded system. In the fol- 
lowing, we will check how the circuit acts on initial states 
|±)^|0);, where the sign can be different on each pair. As we 
will show, all of them give the same output statistics, and thus 
the same holds for their mixture, i.e. the actual initial state. 
These considerations are illustrated in Fig. H] where we take 
the circuits of Fig.[3]and compose them with the decoding and 
encoding circuits (boxed) in order to determine their action on 
the encoded system. 

Firstly, a ZZ{(j)) gate on a pair gives a Z rotation of the 
encoded logical qubit, since the action of ZZ{4)) only depends 
on the parity (Fig.|4^). The action of the second ZZ rotation 
of Fig. [3J) which originally gave an X rotation is shown in 
Fig. ^p. The right hand side is obtained by using CNOT = 
{1®H)CZ{1®H), = 1, the fact that diagonal operators 



2. Toric codes: second scheme 

The second toric-code-like scheme is based on a very differ- 
ent idea. Therefore, observe that the Ky tensor can be written 

as 



K 



v[s\ 



COPY' 



T — H — A[s] — COPY 



(75) 

where COPY is the copy gate |0, 0)(0| + |1, 1)(1|, H is the 
Hadamard gate (both with no physical system associated to 
them), and A the 1-D cluster projector, cf. Eqs. ( [20l ) and ( |2TI ). 
Thus, Ky takes two qubits in correlation space, projects them 
onto the {|0, 0), |1, 1)} subspace, implements the 1-D clus- 
ter map up to a Hadamard, and duplicates the output to two 
qubits. Concatenating these tensors horizontally [this takes 
place in ( |6T1 ) if all A'^'s are measured in Z, and one neglects 
Pauli errors] therefore implements a single logical qubit line, 
encoded in two qubits in correlation space. By removing the 
Hadamard gate from Ky, we obtain a 1-D cluster state en- 
coded in two qubits which is thus capable of implementing 
any one-qubit operation on the logical qubit; in particular, this 
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a) 



b) 



c) 



I'hH — e 



ZZ(<#.) 



[ftA — 



S(4>)- 



:::::::: ^ zz(4,) ^ 




I';, 



zz(0) 



I'), — e 



H — 




ZZ(0) 




ZZ(<#>) 



FIG. 4: Interpretation of the first toric code scfieme in terms of parity encoded qubits. Thie boxed parts of the circuit decode and encode 
the system, a) Z rotations result in Z rotations in the encoded system, b) X rotations result in X rotations in the encoded system, plus Z 
corrections before and after the rotations in case the s qubit below is j— )^ rather than 1+)^. c) Similarly, the coupling circuit Fig.[3]i results in 
a coupling of the encoded logical qubits, up to the same Z correction on the first logical qubit which depends on the s qubit below in exactly 
the same way. Thus, the Z corrections on each qubit cancel out except for the first and the last, which have no effect due to the initialization 
and measurement in the computational basis. 



includes intialization and read-out. We thus define the tensor 



C. Weighted graph states 



COPYt - A[s] - COPY 



(76) 



In this section, we will consider instances of weighted 
graph states iIsI IstIi forming universal resources. To motivate 
the construction, recall that the cluster state can be prepared 
by applying a controlled-phase gate 



Then, the toric code state (|6T] i with Ky replaced by Ky is 
universal for MBQC: Initialization, one-qubit operations, and 
read-out are done exacly as in the 1-D cluster state. The logi- 
cal qubits are decoupled up to Z by-product operators in cor- 
relation space by measuring the Kh tensors in the Z basis. 
The Z by-products in correlation space correspond to Z er- 
rors on the encoded logical qubits and thus can again be dealt 
with as in the cluster. In order to couple two logical qubits, 
we measure a Kh tensor in the Y basis and obtain a ZZ con- 
trolled phase gate in correlation space, which translates to the 
same gate on the logical qubits. Note that this model has the 
additional feature that as as many controlled phases (between 
nearest neighbors) as desired can be implemented simultane- 
ously. 

In the light of the discussion on the initialization of the first 
scheme, one might see similarities between the two schemes, 
since in both cases the information is effectively encoded in 
pairs of qubits. Note however that in the first scheme, the 
information is stored in the parity of the two qubits, and the 
full 4-dimensional space is being used; the reason for this en- 
coding came from the properties of the Kh tensor used as a 
map in horizontal direction. In contrast, the second scheme 
only populates the 2-dimensional even parity subspace, and 
the qubit is rather stored in two copies of the same state; fi- 
nally, the encoding is motivated by the properties of the Ky 
tensor as a map on correlation space in horizontal direction. 



P(0) = |0,0)(0,0| + |0,l)(0,l| + |l,0)(l,0|+e'^|l,l)(l,l|, 

(77) 

with phase = tt between any two nearest neighbors of a 
two-dimensional lattice of qubits initially in the state |+). If 
one wants to physically implement this operation using linear 
optics ll45ll . one encounters the situation that the controlled 
phase gate can be implemented only probabilistically, with the 
probability of success decreasing as (f> increases. It is hence 
natural to ask whether one can build a universal resource using 
gates P{(t>), < < TT, in order to minimize the probability 
of failure* 



1. TranslationaUy invariant weighted grapii states 

Expanding the discussion presented in Ref. iQl], we treat 
the weighted graph state shown in Fig. [5] A tensor network 
representation of these states can be derived along the same 
lines as for the original cluster in Section ITll CI Set \i) := 



* Alternative models with edges resulting from commuting gates with non- 
maximally entangling power can possibly also be constructed by exploiting 
ideas of non-local gates that are implemented with local operations and 
classical communication 14(1. MHi. 
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FIG. 5: Weighted graph state as a universal resource. Solid lines 
correspond to edges that have been entangled using phase gates with 
phase = TT, dotted lines correspond to edges entangled with phase 
gates with (f> = 7r/2. This shows that one can replace some edges 
with weakly entangled bonds. 



2 ^/^(|0) + The relevant tensors are given by 
X 



A[0] 



I ru I ' I lu \ 



A{\\ 



)™l*>iJ->.(l|id(l|.<i(llr 



, (78) 



(79) 



Indices are labeled ru for "right-up" to Id for "left-down". 
The boundary conditions are |0) for the rw, /m, r-directions; 
1+) otherwise. 

We will first describe how to realize isolated evolutions of 
single logical qubits in the sense of Fig. [T] Again the strat- 
egy will be to measure the sites of one horizontal line of the 
lattice in the X-basis and all vertically adjacent systems in 
the Z-basis. The analysis of the situation proceeds in perfect 
analogy to the one given in Section lTlIEl One obtains 



A[Z,+i,u] 



i7S'2^-+^% (80) 



A[Zi+i^d] 



where 



(81) 



and S := diag(l, i) denotes the tt/A gate. 

The operators H and 5* generate the 24-element single qubit 
Clifford group. Following the approach of Section ITVAI we 
take this as the model's by-product group. 

Now choose some phase (p. Re-doing the calculation which 
led to Eq. ( ISOl l. where we now measure in the {|0) ± e"^|l)}- 
basis instead of X on the central node, shows that the evo- 
lution of the correlation space is given by S{(f>), up to by- 
products. In complete analogy to Section HV Al we see that 
the model allows for the realization of arbitrary SU{2) oper- 
ations. 

How to prepare the state of the correlation system for a sin- 
gle horizontal line and how to read read it out has already 



been discussed in Section IIIICI Hence the only piece miss- 
ing for universal quantum computation is a single entangling 
two-qubit gate. 

The schematics for a controlled-Z gate between two hori- 
zontal lines in the lattice are given below. We implicitly as- 
sume that all adjacent sites not shown are measured in the 
Z-basis, 



A[X] 



A[X] 



A[X] 



A\Y] 



A[X] 



A[X] 



(82) 



A[X] 



The measurement scheme realizes a controlled-Z gate, where 
the correlation system of the lower line carries the control 
qubit and the upper line the target qubit. 

In detail one would proceed as follows: first one per- 
forms the X-measurements on the sites shown and the Z- 
measurements on the adjacent ones. If any of these measure- 
ments yields the result "1", we apply a Z-measurement to the 
central site and restart the procedure three sites to the right. 
This approach has been chosen for convenience: it allows us 
to forget about possible phases introduced by other measure- 
ment outcomes. Still, the "correct" result will occur after a 
finite expected number of steps, so the overhead caused due 
to this simplification is only linear It is also not hard to see 
that most other outcomes can be compensated for - so for 
practical purposes the scheme could be vastly optimized. 

Now assume that all measurements yielded "0". Then a Y- 
measurement is performed on the central site, obtaining the 
result y. As we did in Section HV A4| we assume that the 
(lower) control line is in the basis state |c), for c £ {0, 1}. 
The contraction of the lower-most three tensors gives 



JL 



A[X] - A[X] - A[X 



± 



(83) 



5" 



H\c) 



where as before S = S{i) = diag(l,i). We plug this result 
into the A\Y] tensor: 



A[Y] 



(84) 



= \+)lu\+)ru + {-^r^'^{S®S)\+),J+)^ 

Lastly, for x e {0, 1}, 



A[X] — A[X] — A[X] 



= HZ"" 



s-\+) 



5^1 



Hence, the evolution on the upper line is 

H{i + {-ly+yiz) cx Hszy+% 



(85) 



(86) 
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whereas the other one is nothing but the familiar one for a 2-D 
cluster state as in Eqs. (l88l[89] l. with boundary conditions 



\L) = \D) 



\R) = \U) = \l). 



(90) 



The resulting state is hence again a weighted graph state, 
where in one dimension every second edge is replaced by 
an edge prepared using a gate with non-maximal entangling 
power. Then, it is not difficult to see that, again with 

X,Zr,Zu,Zd,Zi e {0,1}, 



FIG. 6: Weighted graph state where the gate is achieved by appropri- 
ately bringing two wires together in a "rerouting process". 



equivalent to Z'^ up to by-products, 
following conclusion: 



We arrive hence at the 



Observation 7 (Non-maximal entangling power). Universal 
resouces may be prepared using commuting gates with non- 
maximal entangling power 



and 



B[X] 



A[Za] 



(91) 



2. Rerouting 

we will consider a second weighted graph state to exem- 
plify yet another novel ingredient that one can make use of 
in measurement-based quantum computation: One can think 
of quantum information being transported in the correlation 
system of some systems on the lattice forming "wires", in a 
way that gates are realized by bringing the "wires" together 
This is an element that is not present in the original one-way 
computer The subsequent example of a resource state has not 
been chosen for its plausibility in the preparation in a physical 
context, but in a way such that this idea of "rerouting quantum 
information" can very transparently be explained, see Fig.|6l 

The resource that we think about is defined by tensors that 
are fully translationally invariant in one dimension, and has 
period two in the orthogonal dimension. 



A — A 



B 



B 



B 



B 



B 



(87) 



A — A 



This is, we have two kinds of tensors: One set is given by 



B[Q] 



I r\ ' I u 



B[\] 



-)r\i)u{W\d 



(88) 



(89) 



A[X] 



(92) 



B[Zd] 



Similarly, we can consider several corner elements in this re- 
source. We obtain 



A[X] - A[Zr] 



i/Z^+^-'S'^", (93) 



B[Zd] 



and similarly 



A[Zi] ^ A[X 



{HSHy''X'"+'', (94) 



B[Zd] 



A\Zu 



T 



B[X] — B[Zr] 



(95) 



A[Zu] 



B[Zi] - B[X] 




(96) 



where we have again made use of the convention that x = 
corresponds to |+) and x ~ 1 to |— ). We need one more 
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ingredient to the scheme, this is 



JL 



and 



B[Zi] - B[0] 



and 



B[Zi] - B[l] 
^ 



)„(0U, (97) 



= \-)r\i)u{l\d, (98) 



A[0] - A[Zr] 





(99) 



A[l] - A[Zr] 





= (-1)"1->«(1|/(1U. (100) 



Putting these ingredients, and following an argument similar 
to the last subsection, we find that up to Clifford group by- 
products, we can transport along the horizontal lines for both 
kinds of local tensors. We can also use the corner pieces to 
reroute as depicted in Fig.|6] and bring routes together forming 
a "gate" imprinted in the lattice, actually, a controlled-S* gate. 

It should be noted that it is not obviously possible to faith- 
fully transport one qubit of information vertically through the 
resource. Loosely speaking, the entanglement between a site 
of type B and the site of type A directly above it is non- 
maximal (this is indicated by dotted lines in Fig. |6]l. Inter- 
estingly, one can still perform a (non-maximally entangUng) 
non-local gate over this connection. 

Observation 8 (Rerouting). Gates in measurement-based 
quantum computation can be achieved by means of appropri- 
ate routing of quantum information in the lattice. 



D. A qubit resource with non-vanisliing correlation functions 

We will very briefly sketch a matrix product state on a 1- 
D chain of qubits, which i) exhibits non-vanishing two-point 
correlation functions, ii) allows for any unitary to be realized 
in its correlation system and iii) can be coupled to a universal 
2-D resource in a way very similar to the AKLT-type exam- 
ple (Section flV At . The discussion will be somewhat superfi- 
cial - however, given the extensive discussion of other models 
above, the reader should have no problems filling in the de- 
tails. 

Choose an integer m > 2 and define 



G := cxp(i7r/mX). 



(101) 



Up to a constant, G is a m-th root of X. The state is defined 
by the following relations: 



|L>=Gt|+), |i?) = |+> 



(103) 



The two-point correlation functions for measurements on 
this state never vanish completely. Indeed, in Appendix 
IVllIAl it will be shown that 



where ^ 2s'm^{TT/m) — 1. 
For X -measurements, we find 



A[X] 



(104) 



(105) 



A[s] 



(102) 



Pursuing the strategy introduced in Section lTV A2[ we set the 
by-product group to B = {Z, G), so the group generated by Z 
and G. One can easily verify that B is indeed a finite group, 
equivalent to the dihedral group of order 2m. 

It is now straight-forward to check that i) measurements in 
the computational basis can be used for preparation and read- 
out (as in Section lTll CI ), ii) general local unitaries can be real- 
ized by means of measurements in the equatorial plane of the 
Bloch sphere (as in Section |TV Alt and iii) a 2-D resource is 
obtainable in a fashion similar to the one presented in Section 
IIV A4I With similar methods, one can also find qubit resource 
states that have a local entropy smaller than unity. 



E. Percolation ideas to make use of imperfect resources 

For completeness, we mention yet another kind of resource: 
This is an imperfect cluster state where some edges are miss- 
ing. Such a setting is clearly relevant in a number of phys- 
ical situations: If the underlying quantum gates building up 
the cluster state are fundamentally probabilistic, such as in 
linear optical architectures, then one very naturally arrives at 
this situation when one aims at minimizing the need for feed- 
forward. A similar situation is encountered in cold atoms in 
optical lattices, when in a Mott state exhibiting hole defects 
some atoms are missing. We do not present details of such 
arguments, which have been considered in Ref. 1431, based 
on ideas of edge percolation and renormalization ||43|l. We 
merely state the result for completeness. Note also that results 
that may be similar to these ones have been announced in Ref. 

We consider the setting where one starts from a 2-D or 3- 
D cubic lattice of size n x n. Two neighboring vertices on 
the lattice are connected with an edge with probability p. The 
stochastic variables deciding whether or not an edge is present 
are assumed to be uncorrelated. If p > p2 — 1/2 holds, then 
it is not difficult to see that one can extract a 2-D renormalized 
lattice of smaller size: This means that one can find a function 
n 1-^ m{n), such that one arrives at a cubic 77i{n) x m(n) ar- 
ray almost certainly as n ^ oo, with the following property: 
Within each of the elements of this array, there is a central site 
that is connected to the central site of the neighboring array. 
Since all the additional sites can be removed by means of Z- 
measurements, we can treat this resource effectively as a 2-D 
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FIG. 7: Cubic lattice of a graph state corresponding to the situation 
where some edges are missing in a cluster state. If the probability p 
of having an edge is sufticiently high the processes independent, then 
a renormalized perfect sublattice can be found almost certainly, giv- 
ing rise to a cluster state of smaller size. If p > P2 = 1/2, where p2 
is the percolation threshold for edge percolation in 2-D cubic lattices, 
then a renormalized lattice can be found almost certainly. Interest- 
ingly, even if 1/2 > p > pz, ps — 0.249 denoting the percolation 
threshold in 3-D, one can almost certainly construct a perfect sublat- 
tice, using an overhead that is arbitrarily close to being quadratic. 

cluster state of dimension 771(71) x m(n), and refer to this as 
a perfect sublattice. This state will not necessarily be exactly 
a cluster state, as it may contain vertices having a vertex de- 
gree of three, but which will nevertheless function as a graph 
state resource just as the cluster state does (for details, see Ref. 
ii^l). Also, n/m{n) is arbitrarily close to being linear in n 
asymptotically. However, an even stronger statement holds: 

Observation 9 (Percolation). Whenever p > p^ = 0.249, /or 
any £ > 0, one can find a function n 1-^ m(n) with the follow- 
ing property: Starting from a sublattice of a 3-D cubic lattice 
of size n X n X 2n/m{n), one can almost certainly prepare 
a perfect sublattice of size 777(77) x m(n). The asymptotic be- 
havior of m can be chosen to satisfy 

n/m{n) = 0{n^). (106) 

That is, with an overhead that is arbitrarily close to the op- 
timal scaUng, one can obtain a perfect resource state out of 
an imperfect one, even if one is merely above the percola- 
tion threshold for a three-dimensional lattice, and not only for 
the two-dimensional lattice, see Fig. [T] The latter argument 
is technically more involved than the former, for details, see 
Ref. ||42|]. This shows, however, with methods unrelated to 
the ones considered primarily in the present work, that also 
random aspects in the resource as such can be dealt with. 



V. ONE-WAY COMPUTATION USING ENCODED 
SYSTEMS 

In the final section of this work, we will show that one can 
find resource states for MBQC that differ substantially from 
the cluster in various entanglement properties. This will be 
done by encoding each system of a resource into several phys- 
ical particles. We will not develop any new computational 
models and make no use of the computational tensor network 
formalism introduced before. The study of encoded resource 



states was initiated in Ref. ID] and later pursued more system- 
atically in Ref. ll25ll . 

More concretely, the following statements will be proved: 

Observation 10 (Resources with weak capabilities for state 
preparation). There exists a family of universal resource states 
such that 

• The local entropy of entanglement is arbitrarily small, 

• The localizable entanglement is arbitrarily small 
and, more strongly, 

• The probability of succeeding in distilling a maximally 
entangled pair out of the resource is arbitrarily small, 
even if one does not a priori fix the two sites between 
which the pair will be established. 

In particular, the resource cannot be used as a state prepara- 
tor. 

We start from a cluster state vector on rz x 77 systems, de- 
noted by I C/„x)i), referred to as logical qubits. As in Ref. [Ql], 
we want to "dilute" the cluster state, i.e. encode it into a larger 
system, by means of invoking the codewords 

|0) := |0)®^ |i) \Wu) (107) 

for some parameter k. The argument relies only on the choice 
of \Wk) as a code word in that we focus on its implications 
on the localizable entanglement, and for that argument, the 
state vector \ Wk) has the desired properties of small local en- 
tropy and permutation invariance. However, for encoded one- 
way computation to be possible, any state vector orthogonal 
to jO)'^'' may be taken, compare also Ref. ll25ll . Every qubit of 
the cluster is subjected to the encoding operation 

|0)(0| + |i)(l| (108) 

yielding the diluted cluster \Dn.k)- A set of physical qubits 
corresponding to one cluster bit will be called a block. As be- 
fore, by a local measurement scheme we mean a sequence of 
adaptive local projective measurements, local to the physical 
systems. 

Let us first show again in more detail that such an encod- 
ing constitutes no obstacle to universal quantum computation. 
Each of the code words is orthogonal, and for computation 
to be possible, we need to do local dichotomic measurements 
in the logical space. By Ref. fiill, any two pure orthogonal 
multi-partite states on k qubits can be deterministically dis- 
tinguished using LOCC. By making use of the construction 
of Ref. ll44ll . this can be done by an appropriate ordered se- 
quence of adapted projective measurements tti ® • • • (X) Tr^ on 
the sites of each codeword, giving rise to an arbitrary projec- 
tive dichotomic measurement with Kraus operators 

A2 := = 1 - IV) (^1 (109) 

in the logical space, = a|0)+/3|l) and \i>^) = -/3*|0) + 
q:*|1). Hence, one can translate any single-site measurement 
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on a cluster state into an LOCC protocol for the encoded clus- 
ter. This shows that j^P) is universal for deterministic MBC. 
This is the argument of Ref. (HI] (see also Ref. ll25ll for a 
more detailed and extensive discussion on one-way comput- 
ing based on encoded systems). 

In the following we are going to show in more detail that 
despite this property, we are heavily restricted to use this re- 
source to prepare states with a significant amount of entangle- 
ment between two constituents. In fact, we can not even dis- 
till a perfect maximally entangled qubit pair beyond any given 
probability of success. This means that these states are univer- 
sal resources, but on the level of physical systems utterly use- 
less for state preparation. The given resource is, needless to 
say, not meant as a particularly feasible resource. Instead, we 
aim at highlighting to what extent as such the entanglement 
properties can be relaxed, giving a guideline to more general 
settings. 

Note first that the localizable entanglement El in these 
resources can easily be shown to be arbitrarily small: The 
entropy for a measurement in the computational basis reads 
Hb{3/{4k + 2)), where Hb : [0, 1] [0, 1] is the standard 
binary entropy function. Using the concavity of the entropy 
function, we find 

EL{\V„^k){Vn.k\) < Hbi3/{4k + 2)), (110) 

such that \imk-foo ELdVn^k) {T^7i,k\) = 0. This means that 
for two fixed sites, the rate at which one can distill maximally 
entangled pairs by performing measurements on the remain- 
ing systems is arbitrarily small. 

This can be seen as follows. We will aim at preparing a 
maximally entangled state between any two constituents of 
two different blocks. It is easy to see that within the same 
block, the probability of success can be made arbitrarily small. 
We hence look at a LOCC distillation scheme, a measurement- 
based scheme, taking the input p and producing outputs 

p^KjpK] (111) 

with probability pj = ti{KjpKj), j = 1, . . . , J. This cor- 
responds to a LOCC procedure, where each of the measure- 
ments may depend on all outcomes of the previous local mea- 
surements. Let us assume that outcomes labeled 1, . . . , S for 
some iS* < J are successful in distilling a maximally entangled 
state. 

We start by exploiting the permutation symmetry of the 
code words. Choose a block i of |I?„,fe). Assume there ex- 
ists a measurement-based scheme with the property that with 
probability p, the scheme will leave at least one system of 
block i in a state of maximal local entropy. Then there exists 
a scheme such that with probability p, the scheme will leave 
the first system of block i in a state of maximal local entropy. 
At some point of time the scheme is going to perform the first 
measurement on the i-th block. Because of permutation in- 
variance, we may assume that it does so on the fc-th system of 
the block. The remaining state is still invariant under permu- 
tations of the first fc — 1 systems. Hence there is no loss of 
generality in assuming that the next measurement on the i-th 



block will be performed on the k — 1-st system. If the local 
entropy of any of the unmeasured systems is now maximal, 
then the same will be true for the first one - once again, by 
permutation invariance. 

Also, it is easy to see that the probability p that a 
measurement-based scheme will leave any system of block i 
in a locally maximally mixed state is bounded from above by 

p<2/k. (112) 

Let pi be the initial probability of obtaining the outcome |1) 
for a Z measurement on this qubit, pi = Clearly, 

PiKl/k. (113) 

We consider now a local scheme potentially acting on all 
qubits except this distinguished one, with branches labeled 
i = 1 , . . . , J, aiming at preparing this qubit in a maximally 
mixed state. Let ps be the probability of the qubit ending up 
in a locally maximally mixed state. In case of success, so in 
case of the preparation of a locally maximally entangled state, 
we have that pi (s) = 1/2, in case of failure pi (/) > 0. Com- 
bining these inequalities, we get 

l/k > = p,pi (s) + (1 - p,)pi (/) =pj2. (1 14) 

We can hence show that there exists a family of universal re- 
source states such that the probability that a local measure- 
ment scheme can prepare a maximally entangled qubit pair 
(up to l.u. equivalence) out of any element of that family is 
strictly smaller than e > 0. 

Let Pi be the probability that a site of block i will end up as a 
part of a maximally entangled pair. This means that when we 
fix the procedure, and label as before all sequences of mea- 
surement outcomes with j ~ 1, J, one does not perform 
measurements on all constituents. Let / denote the index set 
labeling the cases where somewhere on the lattice a maxi- 
mally entangled pair appears, so the probability p for this to 
happen is bounded from above by 

iG/ 

According to the above bound, Pi < 2 /fc, giving a strict upper 
bound of p < 2n^ /k for the overall probability of success. 
The family 

1*0 := |25n.fc(n)>' (116) 

for k{n) 2e^^n'^ is clearly universal, involves only a linear 
overhead as compared to the original cluster state and satisfies 
the assumptions advertised above. 



VI. CONCLUSIONS 

In this work, we have shown how to construct a plethora of 
novel models for measurement-based quantum computation. 
Our methods were taken from many-body theory. The new 
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models for quantum computation follow the paradigm of lo- 
cally measuring single sites - and hence abandoning any need 
for unitary control during the computation. Other than that, 
however, they can be quite different from the one-way model. 
We have found models where the randomness is compensated 
in a novel manner, the length of the computation can be ran- 
dom, gates are performed by routing flows of quantum infor- 
mation towards one another, and logical information may be 
encoded in many correlation systems at the same time. What 
is more, the resource states can in fact be radically different 
from the cluster states, in that they may display correlations 
as typical in ground states, can be weakly entangled. A num- 
ber of properties of resource states that we found reasonable to 
assume to be necessary for a state to form a universal resource 
could be eventually relaxed. So after all, it seems that much 
less is needed for measurement-based quantum computation 
than one could reasonably have anticipated. This new degree 
of flexibility may well pave the way towards tailoring compu- 
tational model towards many-body states that are particularly 
feasible to prepare, rather than trying to experimentally realize 
a specific model. 



so for s e {0, 1}: 

+ (119) 

where := 1,1 := and p := |(0|G|0)p = sin^(7r/TO). In 
other words: when acting on the computational basis, $ im- 
plements a simple two-state Markov process, which remains 
in the same state with probability p and switches its state with 
probability (1— p). Now, {ZiZi+k) equals +2 if an even num- 
ber of state changes occurred and —2 if that number is odd. So 
for the expectation value we find 

= 2Y,(%'^-\l-pf{~lf (120) 
(=0 ^ ^ 

= 2(2p- 1)'' = 2(2sin2(7r/m) - 1)*=. 
B. Hamiltonian of the AKLT-type state 
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VIII. APPENDIX 

A. Computing correlations functions 

What is the value of the two-point correlation function 
{ZiZi^k) — {Zj) {Zi+k)'^ In this work, we have only intro- 
duced the behavior of the correlation system when subject to a 
local measurement of a rank-one observable. However, in or- 
der to evaluate the correlation function, we need "measure the 
identity" on the intermediate systems or, equivalently, trace 
them out. Without going into the general theory ||26ll , we just 
state that tracing out a system will cause the completely posi- 
tive map 

^ : p^^A[i]pA[i]'^ (117) 

i 

to act on the correlation system. 

For the cluster state, using the fact that the bases {|0), |1)} 
and {|+), |— )} are unbiased, we can easily show that $^ is 
the completely depolarizing channel, sending any p to 2~^1. 
This causes any correlation function to vanish for k > 2. How 
does the situation look like for the state vector defined by Eq. 
( fT02l i? We compute: 

$:pK^ ^ tr(pG|s)(s|Gt)|0)(0|, (118) 

s=0,l 



In Section ITV Al we discussed an AKLT-type matrix prod- 
uct state. It was claimed that the state constitutes the unique 
ground-state of a spin-1 nearest neighbor frustration free 
gapped Hamiltonian. It must be noted that in this work, we 
have not introduced the technical tools needed to cope with 
boundary effects at the end of the chain. There are at least 
three ways to make the above statement rigorous: a) treat the 
statement as being valid asymptotically in the limit of large 
chains, b) work directly with infinite-volume states ll26ll . or 
c) look at sufficiently large rings with periodic boundary con- 
ditions iZTll . Once one chooses one of the options outlined 
above, the proof of this fact proceeds along the same lines as 
the one of the original AKLT state, as presented in Example 7 
of Ref. It26,l (see also Ref. li27ll ). Indeed, using the notions of 
Refs. iloL I27I1 one verifies that 

La : 6(C2) ^ C2, (121) 
B ^ EL»2=i HBA[H]AM\iui2) (122) 

is injective. Further, if := range it is checked by direct 
computation that dim(CJ2 (8) 1 n 1 (g) f/2 ) = dim 02- All claims 
follow as detailed in Refs. l|26i|27ll . 

In particular, let /i be a positive operator supported on the 
vector space spanned by: 

{|1,1),|2,2),-(1/4)|0,0) + |1,2) + |2,1), (123) 
-(l/^/8)|0,0) + |0,2) + |2,0), 
-(1/V8)|0,0) + |0,1) + |1,0)}. 

Set H where translates its argument i sites 

along the chain. Then iJ is a non-degenerate, gapped, frustra- 
tion free, nearest neighbor Hamiltonian (called parent Hamil- 
tonian in Ref. ll27ll ). whose energy is minimized by the state 
at hand. 
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